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A Zeeman field affects the metallic heavy fermion ground state in two ways: (i) it splits the 
spin-degerenate conduction sea, leaving spin up and spin down Fermi surfaces with different band 
curvature; (ii) it competes with the Kondo effect and thus suppresses the mass enhancement. Taking 
these two effects into account, we compute the quasiparticle effective mass as a function of applied 
field strength within hybridization mean field theory. We also derive an expression for the optical 
conductivity, which is relevant to infrared spectroscopy measurements. 



I. INTRODUCTION 

Heavy fermion materials^ are metallic alloys of ac- 
tinide or rare-earth elements (typically U or Ce) . At high 
temperature, their chemically active valence electrons are 
confined in localized / orbitals. These /-electrons consti- 
tute a dense lattice of localized spins embedded in — and 
only weakly interacting with — the ordinary conduction 
sea. At low temperature, the /-electron moments be- 
come strongly coupled to the conduction electrons and, 
indirectly, to each other. Below the characteristic Kondo 
temperature, a complicated many-body state emerges in 
which the local moments arc screened (provided that the 
Kondo physics overwhelms the competing RKKY inter- 
action^). This state exhibits unconventional metallic be- 
haviour with an effective mass tens or hundreds of times 
larger than that of a bare band electron. 

The key detail is that the broad band of conduction 
electrons is intersected by a nearly dispersionless, highly 
correlated, /-electron band. The heavy fermion ground 
state can be understood to be a "nearly-broken symme- 
try" state 4 in which the overlap between excitations in 
the two bands becomes macroscopically important (to or- 
der l/N, where N — ► oo represents the limit of large or- 
bital degeneracy 4 *^) . This motivates treating the Kondo 
physics, in the guise of an interband hybridization, at the 
mean field level. 

When the hybridization order parameter condenses, 
the /-electrons are incorporated into an enlarged Fermi 
sea of composite quasiparticles. Mixing between the local 
and itinerant degrees of freedom causes the conduction 
band to break into upper and lower quasiparticle bands, 
leaving a region of shallow dispersion (i.e., large effective 
mass) near the hybridization gap edge. 

The hybridization picture has been verified in real ma- 
terials by a variety of experimental methods. Infrared 
(IR) spectroscopy in particular has proved to be an im- 
portant experimental probe: IR studies of YbFe4Sbi2 
and CeRu4Sbi2 provided the first direct observation of 
the hybridization gap in heavy metals p£ In addition, the 
complex optical conductivity and the dielectric function, 
obtained from reflectance measurements by Kramers- 
Kronig, can be used to infer the effective mass. Such an 
analysis is possible because the electrodynamic response 
of the system depends on the mass enhancement factor 



in a known way&2ii£ 

In a recent experiment, Dordevic et al. carried out a de- 
tailed IR study of the heavy fermion material CeRu4Sbi2 
in high magnetic fieldsiii A complete theory of the elec- 
trodynamic response was lacking in this case. Here, the 
applied field breaks the spin degeneracy of the Fermi sea, 
with the result that spin up and spin down quasiparticles 
have different effective masses. The reflectance, however, 
is insensitive to the quasiparticle spin. Thus, the mea- 
sured optical conductivity, insofar as it depends on the 
(now spin-dependent) mass enhancement factor, repre- 
sents an average over the two spin projections. 

In this paper, we compute the mean field heavy fermion 
ground state in the presence of a Zeeman field (under the 
assumption that the field does not induce a magnetic in- 
stabilitjii2ii£). We find that the effective mass is modified 
to linear order because of Fermi surface splitting and to 
quadratic order because of magnetic suppression of the 
Kondo singlet amplitude. As it turns out, the optical 
conductivity averages over the effective mass in such a 
way that corrections first appear at second order and 
with a negative coefficient; the characteristic field scale 
is on the order of the Kondo energy (in magnetic units). 
Such behaviour is observed experimentally in Ref. Illl 

The paper is organized as follows. The hybridization 
mean field equations are derived in Secs.lTTland lllll Their 
solution in zero and nonzero field is presented in Sec. II VI 
followed by a discussion of the optical conductivity in 
Sec. El 

II. MEAN FIELD FORMALISM 

A. Hybridization Field 

The Kondo Lattice Model (KLM) is thought to provide 
an approximate description of heavy fermion materials. 
(The Periodic Anderson Model in the limit of large onsite 
repulsion is an equivalent starting point.—) It describes 
a band of free conduction electrons moving in a periodic 
array of magnetic impurities S r . The only interaction is 
the onsite Heisenberg exchange between the electron spin 
density and the magnetic moment of the impurities. 

Hklm = -^2t r - r >c\,c r > + J~Y^s r ■ S r . (1) 

rr' r 



2 



Here, t r = t*_ r and J are the hopping and exchange 
integrals, respectively, and (c) is the creation (annihi- 
lation) operator for the conduction electrons; s = ^c^crc 
describes their local moment. The impurity spins are 
S = 1/2. The model can be extended to include a Zee- 
man term, which consists of a magnetic field B coupled 
to the total magnetic moment at each site: 



H = H) 



KLM 



B 



(s r -\- S r 



A useful formal trick is to represent the local spins 
in terms of fermions, which puts the local and itinerant 
degrees of freedom on the same footing. Specifically, we 
take S — \f*trf subject to the constraint f>f = 1. Here, 
f' (/) is the creation (annihilation) operator of a ficti- 
tious, dispersionless / band. The constraint suppresses 
all /-charge fluctuations and has the effect of projecting 
out the singlet states that are not a part of the Hilbert 
space of the original SU(2) spin. We shall assume that 
it is sufficient to enforce the constraint on average and 
require only that (Pf) = 1. 

The operators x M = ^=/ 'a^c, defined in terms of the 

unit and Pauli matrices = (1, <t), describe the spin 
degrees of freedom at each site. They can be used to 
express the exchange term in the Hamiltonian as 



-!x ot x° 



4 X 



(3) 



X° and x ac t m the singlet and triplet channels, respec- 
tively. In the heavy fermion state, the singlet amplitude 
condenses (x°) ^ 0, leading to hybridization of the c 
and / bands. Neglecting second order fluctuations about 



(x°> 



gives 



X 0t X° 



<x }*x° + x ot <x M<x )| 



(4) 



The resulting mean field Hamiltonian is 

HmF = - ^ 4 [*r-r' + 8rr> (Mc + B ■ cr/2 
rr' 

-E#0/+-B-<r/2)/r 

r 



(5) 



We have expressed the Hamiltonian in terms of a hy- 
bridization field 



(6) 



having units of energy. We have also included two chem- 
ical potentials fi c and [it, which couple to the c- and 
/-electron densities. These allow us to control the con- 
duction band filling and the /-level occupancy. 



B. Quasiparticle Dispersion Relation 

The minimum-energy hybridization field configuration 
is translationally invariant. Thus, we can write V r = V = 
\V\e l6 . Moreover, although the hybridization (pc) is in 
general complex, we may safely take V to be real and 
positive. There is a U(l) gauge freedom associated with 
the invariance of S — hptrf under the phase rotation 



(2) / e^f. Accordingly, by fixing 



—tf, we can gauge 
away the phase of V. 

Equation (|SJ| has a particularly simple wavevector rep- 
resentation, 



H MF — 



i.- 



(4 fL) M ks (£) + 



3J 



(7) 



In the summation, fe ranges over all wavevectors in the 
Brillouin zone and s over the two fermion spin projec- 
tions. N denotes the number of lattice sites. The coeffi- 
cient matrix 



M fc 



sB/2 



-V 



-V -Hf-sB/2) 



is a function of the free conduction-electron dispersion 



-k-r 



t r , the physical parameters J and B, 



£fc : 

and the mean field parameters fi c , [if, and V. By writing 
B ■ cr ss i = Berks' = sBS ss i , we have chosen to direct the 
applied magnetic field along the 3 axis. 
It is convenient to introduce the quantities 



2/x = fi c + [if and b = [i c — [if, 



(9) 



which serve as an alternate set of Lagrange multipliers 
(rotated 45 degrees with respect to the original set). As 
we shall see, [i is the energy required to remove a dressed 
quasiparticle from the top of the Fermi sea, whereas 6/2 
is the energy to remove a bare conduction electron. We 
can think of [i as the chemical potential of the fully in- 
teracting system. The parameter b controls how many 
electrons are available to compensate each local spin. 

Let us shift the diagonal entries in Eq. © so as to void 
the lower right entry. The resulting matrix is 



M ks + ([if + sB /2)l 



e k -b -V 
-V 



which has eigenvalues 



r- 1 - 
ii- — — 



e k - b + ny/(e k - b) 2 + AV 2 



(10) 



(11) 



labelled by n = ±. In order to diagonalize Eq. (JHJ, we 
construct a unitary transformation U k , whose columns 
are populated with the normalized eigenvectors of M ks : 



U 



JJC+ JJC- 

u*+ uf- 



-r 



J(i+) 2 +v 2 J(i-) 2 +v 2 

V V 
^(I+) 2 + V 2 y/(I-) 2 +V 2 



(12) 
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Note that all the k dependence of Uk comes from IP and 
that both and Uk are independent of the magnetic 
field. 

One can easily verify that Eq. I|12|) is unitary and that 
it diagonalizes the coefficient matrix: (UlUk) nn — S nn 
and {UlM ks Uk) nn ' = (E% 8 - rfS™' , where 



Wavevector Sum DOS Integral 



l r 



e k - sB + ny/(e k -b)* + W 2 



(13) 



When 7^0, the c- and /-electrons admix to produce the 
quasiparticles of the interacting system. Their dispersion 
is given by Eq. (|13|l . Note that hybridization splits the 
band structure into two disjoint pieces. The quantum 
number n labels quasiparticles in the upper (n = +) and 
lower (n = — ) bands. 

In a Zeeman field, the Fermi surface is split into two 
distinct Fermi surfaces — one for spin up and one for spin 
down quasiparticles. For quasiparticles with spin pro- 
jection s, the Fermi surface is defined by E£ s = Li or, 
equivalently, I£ = fif s = /Sf + sB/2. 



C. Fixing the Mean Field Parameters 

The physics of the hybridization mean field theory de- 
pends on how the parameters fi c , /if, and V vary as a 
function of the Kondo coupling and the magnetic field. 
The optimal value of the hybridization strength is the 
one that minimizes the free energy density 



[-> 

Cfc - 
n 



— V 2 /cu 
{lo 2 - V 2 )/u 
sgn(w) 
/ duD{io) 



TABLE I: A wavevector sum is transformed to a density 
of states integral by making the substitutions listed in this 
conversion chart. In addition E% a = fi corresponds to uj — 



D. Thermodynamic Limit 

In pursuit of a solution to Eqs. (|16[) . it is helpful to 
eliminate the k summations in favour of energy integrals 
weighted by the density of states (DOS), 



kn 



(17) 



In the thermodynamic limit (N ~ > oo), the set of Ij£ val- 
ues is dense, and Eq. (|17fl is a smooth function of w. [As 
a convenience, we have defined the DOS as the spectrum 
of 1% rather than E^ s , which makes the function inde- 

The true DOS is offset 
= fif s - fx.] Applying the 



pendent of the magnetic field, 
from Eq. (HJ by -6/2 + sB/2 
delta function identity 



T = 



8V 2 
3J 



1 

Jn 



5^1n[l + e-^.-") 

hsn 



(14) 



The Lagrange multipliers are chosen to satisfy (c^c) = 
—dT/djic = n c , where n c is the conduction band filling, 
and (f'f) = —dJ^/dfif = 1. Expressing these conditions 
in terms of the rotated Lagrange multipliers of Eq. JjJJ , 
we find that 



-2^=n c -l, and - = 0. (15) 



Performing the free energy differentiations yields 

^£/(^-M) = nc + l (16a) 



ksn 



N ^ 



3 J 



n(e k - b)f(El s - g) 



= n r -l 



= 0. 



(16b) 
(16c) 



Note that Eq. (|16afl defines the Luttinger volume as n — 
n c + 1, with both the c- and /-electrons counted in an 
enlarged Fermi sea. Equation l|16cfl is the gap equation 
that determines V. 



n=±l 



£fe-6-2w + W(£ fe -6) +W 2 



= 2 n 



0J Z 



£k 



V 2 

( 1 - -J 



(18) 



to Eq. (fT7|> . we can show that D(ui) ~ 1 + V 2 /uj 2 . The 
correspondence between quantities in the wavevector and 
energy representations is summarized in Tabic [I] 

The most important feature of the DOS is that it de- 
velops a band gap as V increases from zero. The DOS of 
the noninteracting conduction electrons can be written 
as the product of a line-shape function g(u>) and a heav- 
iside function, which ensures that the density of states 
vanishes outside the band: 

D oM = jf E ^ - «*) = 9(u)6(W 2 - W). (19) 



In the interacting system, the conduction-electron DOS 
has the same basic form, 

^» = ^EW^-tf) = fl(AM)0(W 2 -4AM 2 ), 



ki i 



but its energy scale is renormalized by the function 
X(u>) = u)(l — V 2 /us 2 ) + 6, which comes from the delta 
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FIG. 1: The total density of states for the noninteract- 
ing system Do(ij) consists of an /-level delta function peak 
superimposed on a conduction band of width W. When 
V / 0, hybridization breaks the spectral weight into lower 
(u>\ < W < U02) and upper (103 < ui < U4) bands. 



function on the right-hand side of Eq. I|18(l . As a result, 
the argument of the heaviside function in Eq. (|2(J[1 is a 
fourth degree polynomial whose roots 



LOl = - 


vV/ 4 4- 


-b/2) 2 - 


W 2 


- W/4- 


b/2 


(21a) 


0J 2 = - 


vV/4- 


-b/2) 2 - 


W 2 


f W/4- 


b/2 


(21b) 


^3 = + 


sJ{W/i-\ 


-b/2) 2 - 


W 2 


- W/A- 


b/2 


(21c) 


CJ4 = + 


\J(W/4- 


■ b/2) 2 - 


W 2 


FW/4- 


b/2 


(21d) 



W/2- 



W/A- 



- 



-W/4- 



-W/2- 
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|2A K 
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delineate the band edges: 



FIG. 2: The hybridized band structure is plotted in zero 
field for conduction band filling n c — 3/4 and Kondo energy 
Ak = W/25. The heavy metal has a partially filled lower 
band with quasiparticle states occupied up to the Fermi 
wavevector A;f. Several important energy scales are marked 
with (red) dotted lines. 



9{W 2 - 4A(w) 2 ) = ^(-1) i+1 0(cj - ut). 



(22) 



III. HEAVY FERMION METAL 



Spectral weight exists only at energies in a lower band 
from uj\ to L02 and in an upper band from L03 to L04. The 
two bands are separated by a gap of width 2 A = u)% — 102- 
See Fig. □ 

A similar analysis shows that the DOS for the /- 
electrons differs from D c (lu) by a factor of V 2 /lu 2 . Hence, 
D(uj) = D c {uj) + D f (uj), the total DOS, is equal to 

D{u) = (l + K)g{X{u)) £(-l) l+1 ^ - <*)■ (23) 

Using the conversion chart in Table we can re-express 
Eq. JT3J as 



T= S I1-\Y f duJ D{uj) 1J1 + e -K^fs) 
3 J P ^ J L 



(24) 



Hence, the constituent equations of the mean field theory 
[viz., Eqs. 1)16(1 ] can be written compactly as 











Df(cj) 


hi 


[:! 


> . (25) 











A. Characterizing the Ground State 

The key feature of the heavy fermion state is the hy- 
bridization gap. It generates a region of very shallow 
dispersion near the gap edge, which is responsible for the 
large effective mass of the quasiparticles. For concrete- 
ness, let us suppose that the lower band is filled to some 
point below the hybridization gap. Then the relevant dis- 
persion relation is that of the lower band. This situation 
is depicted in Fig. [5] The heavy fermion state is metallic 
and possesses a well-defined Fermi surface given by the 
set of fc points satisfying = fi. 

The effective mass to* of the quasiparticle excitations 
is a function of the band curvature. It is related to the 
nonintcracting band mass by the variation m*5Ej^ s = 
m*5Ifc = mSek, averaged over all points on the Fermi 
surface. Hence, the mass enhancement factor is given by 
(m*/m) = {dlji/dek^ 1 and, via Tabled 



m*\ d ( l 



1 v 2 
= l + — 



(26) 



The direct energy splitting between the bands can be 
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written as 



^ks ^ks <2 



(27) 



Since optical experiments probe quasiparticles at the top 
of the Fermi sea by promoting them to the upper band 
with negligible momentum transfer and without inducing 
spin flips (dk = 0, s fixed), the optical gap is determined 
by evaluating Eq. I)27|) at the Fermi level. Again, via 
Table H 



/ . \ 1 , ,/ V 2 
(A opt j s = - sgn(w) I uj + — 



(28) 



Note that Eqs. (|26fl and (|28|l are equal up to an overall 
factor proportional to the /-level energy: 



iopt 



m* /m 



(29) 



B. Mean Field Equations in Detail 

Let us represent the bare conduction band filling by 
n c = 1 — x. (The volume enclosed by the enlarged Fermi 
surface is n c + rif = 2 — x.) We can think of x as the den- 
sity of holes in the lower band doping the system away 
from the half-filled Kondo insulator state. Further, since 
Van Hove singularities do not play an important role 
here, let us assume that the density of levels in Eq. 
is fiat and replace the line-shape function by its average 
value g -> l/W. The total DOS is then 



D(u>) 



' / i + ^)^(-ir +i ^- w . i ). (so) 



W 



Expanding Eqs. in V 2 /(W±b) 2 < 1, we find that 
the bottom of the band is given by 



W 2 



(W 2 f 



UJl 2 2(W + 2b) ' 2(W + 2b) 3 

and the hybridization gap edges by 

W 2 (W 2 ) 2 



UJ 2 



U) 3 = 



2(W~2b) 2(W-2b) 3, 

W 2 (W 2 ) 2 
2(W + 2b) ~ 2(W + 2bf 



(31) 

(32a) 
(32b) 



Defining the gap width 2 A = C1/3 — ci>2, we can invert 
Eqs. (|32|l to solve for the hybridization strength: 



W 2 = 2A 



(W + 2 A) -(W- 6A) 



(33) 



Since b controls the difference between the c- and /- 
electron occupation, it is proportional to the number of 
holes in the lower band. We can express b = —Wx/2 + e 



as the noninteracting (J = 0) result plus a correction 
e ~ A. Equations and (gSJ become 

4Ae (2A) 2 

-2wi = W(l-x)+2e+2A(l+x)- — -^-2x, (34) 



-2lo 2 



2A 1 



4Ae 



(2A) 5 
W 



2x, (35) 



and 



AV 2 = 2A W(l - x 2 ) + 4.xe + 2A(l + 3.x 2 



(36) 



From here, it is straightforward to sketch out how the 
solution to the mean field equations is obtained. The 
zero-temperature c- and /-electron occupation are com- 
puted by integrating D c (ui) and (ui) from u\ up to 
fif s . So long as /x/ s < 102 , we can write 



(etc) 



W 



E 



du) 



W 



(Hf-wt) (37) 



and 



W 2uj 1 {^ 2 -B 2 ) 



(38) 

The first condition in Eq. I|25|) . (e'e) = 1 — x, implies that 
2[if — 2u>i = W(l — x), which defines — fif ~ e + A(l + x) 
via Eq. I|34|l . The second condition, (f'f) = 1, fixes the 
value of e in terms of the variable A (the only remaining 
unknown) and the constants x and W. Finally, the third 
condition, 



1 = 



3 J 



3J 

aw 



in 



An 2 - B 2 



W(l-x)-2fi f 



(39) 

closes the system of equations. In the next section, we 
carry out these steps explicitly for the B = and B ^ 
cases. 



IV. MEAN FIELD SOLUTION 



A. In Zero Applied Field 

When B = 0, the requirement that Eq. (|3*8f equal 
unity reduces to 

AV 2 (l-x) =4m/^i = -2wi [-2u) X - W(l - x)] . (40) 

Substitution of Eqs. I|34|) and then allows us to solve 
for e as a function of the hybridization gap: 

(2A) 2 

2e = -2Az(l + x) + ^—^-x 2 (l - x). (41) 
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Several results follow immediately. The energies of the bottom and top of the lower quasiparticle band are 

2A(W)-( 2A ) 2 
(2A) 2 



-2wi 
-2uj 2 



The Lagrange multipliers arc 



/'<-: 



W{l-x) 
2A(l-a;' 



W 



-x(l 



W 



'(1 



(2A) 2 



W 



(1 



2Aa 



(2A) 2 
W 



(42) 
(43) 

(44) 
(45) 



or, alternatively, 



= - A- 



VF + 2A(l 



(2A) 5 



ir 



(1 - a;) (2 - x) 



W^ + 2A(l + a;) 



(2A) S 
VF 



x(l — an 



(46) 
(47) 



The hybridization energy is 



AV 2 = 2A(1- 



W(l 



2A(l + a: + 2.T 2 ) 



r 



(48) 



As a consistency check, we verify that the assumptions 
made during the derivation hold true: the /-level chem- 
ical potential, /if = — A(l — x 2 ), does indeed sit below 
the top of the lower band, 0J2 = — A(l — x); it is also true 
that b ~ x and e ~ A. 

The value of A is obtained by substituting Eq. 145(1 
and B = into Eq. 



3J 



In 



2A{l-x 2 



AW \W(l- x) +2A(1 



(49) 



The solution of this gap equation is the Kondo energy, 



2A K = 



1 + x (l + x) 3 



(50) 



For realistic values of the physical parameters, the bare 
exchange coupling is smaller than the bandwidth. Thus 
a = e - 4M 7 3J < 1, which implies that 2A K /W, 



-2nf _ a(l - x) 



W 



1-a 



and 



AV 2 
W 2 



= a(l + 2a)(l-a;) (51) 



are small parameters. The expansion leading to Eqs. 1(31(1 
and 1(32(1 . however, requires the stronger condition that 



AV 2 



(W — 2b) 2 



1 -x 



(52) 



be small. The hybridization mean field theory is not ap- 
propriate when n c — 1 — x < a, the so-called "exhaustion 
limit" of Nozieres^ 



Most important, small a guarantees that the mass en- 
hancement factor is a large number: 



m 
m 



1 



a(l — x) 



1. 



(53) 



We emphasize again that the enhancement is a conse- 
quence of the small hybridization gap and the very shal- 
low quasiparticle dispersion at the top of the Fermi sea, 
as depicted in Fig.[2J The energy scale of the optical gap, 
also shown in the figure, is the conduction bandwidth, a 
appears only as a subleading contribution: 



a _ W 

A op t — — 



1 



(2-x) 



(54) 



B. In Nonzero Applied Field 

When B =/= 0, the quasiparticles acquire a net magne- 
tization, M = i((ct(T 3 c) + (/V 3 /)). Since 



(cVc) 



1 

W 



and 



AV 2 



B 
W 



B 



W Afi 2 



B 2 



(55) 



(56) 
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FIG. 3: The plot shows the band structure of Fig. Sunder the 
influence of an applied Zeeman field. The quasiparticle bands 
E^ s are spin split, and the Kondo energy Ak is somewhat 
reduced from its zero-field value. The inset (top-left) pro- 
vides a magnified view of the (dotted) region near the lower 
hybridization gap edge. The difference in band curvature at 
&Ff and fcp| is clearly visible. 



the spin susceptibility \ = dM/dB\B=o is Pauli-like, but 
enhanced by a factor oil/ a with respect to that of the 
bare conduction electrons [cf. Eq. (|26|l ]: 



X 



2W 



Xc 



(57) 



In terms of the solution outlined in the previous sec- 
tion, the effect of the applied field appears as a modifica- 
tion to Eq. (|4*T|) . Imposing (c^c) — 1 — x and {ft f) = 1 
by way of Eqs. (JSZJl and (jSHl yields 



2e = ^/A 2 (l-x 2 ) 2 + B 2 - A(l 



(2A) ; 
W 



l (l-x) + 



B 2 



■Ax + l 



W {l + x) 2 {l-x)' 



(58) 



This result is correct to up to terms of order B 4 /W 3 , 
which are negligible in the regime where the applied 
field is comparable to the Kondo energy (B < A <C 
W). To second order in B/A, Eq. ifSSI) behaves as 
2e = (2e) B =o + B 2 /2A(1 - x 2 ), and thus induces field- 
dependent corrections to Eqs. I|42 |) -i|18 )l in the obvious 
way. In particular, we have 



-2fi f = 2A(1 - x 2 ) + 



B 2 



2A{l-x 2 ) 



(59) 



Once again, the system of equations is closed by ap- 
pealing to Eq. (|39|1 . which in nonzero field admits the 
series solution 



-2 H 



aW{l 



B 2 



(1 - a 2 )B 4 



1-a 2aW(l - x) %aW[l - x) 
To leading order in a, comparison with Eq. I|59[) gives 



(60) 



2A K = 



aW 
1 + x 



1 - 



1 I B 



2\B t 



(61) 



Here, we have introduced the characteristic field strength 
Bo = /i/|s = o = ccW(l — x). Note that the enhancement 
of the /-level chemical potential is accompanied by a sup- 
pression of the Kondo energy. This is a consequence of 
the Zeeman field's favouring triplet over Kondo singlet 
pairing. 

The corresponding results for the effective mass and 
optical gap are 



rn 
rn 



rn 
rn 



B\ 3 / B 
~ S \ Bo J + 2 I B~ 



and 



(A„ p t) a = (A opt ) 



B 
Bo~ 



XQ 

~2~ 



B 



(62) 



(63) 



In Eq. (|62|l . the linear correction arises from the (if + 
sB/2 denominator in Eq. (|26|l and can be directly at- 
tributed to the spin splitting of the Fermi sea. The 
quadratic correction incorporates both Fermi-sea effects 
appearing at second order and Kondo energy suppression 
appearing indirectly though the V 2 numerator. Equa- 
tions (|f]2l and $53$ satisfy the ratio rule given in Eq. 



V. OPTICAL CONDUCTIVITY 

The total current density in the presence of a vector 
potential A follows from the Kubo formula in the usual 
way: 



e 2 1 v ^ 



ksn 



dE, 



ks 



to 



v k -q + iT 



2 gap 



-A^q,u). (64) 



Here, v% = dek/dk a is the velocity of the bare band 
electrons and 1/r is the quasiparticle lifetime. The cur- 
rent density is related to the conductivity by Ohm's law, 
J = uE = aiuiA/c. Hence, 



a a0 {q 



0,,) = ^= lim^- f^'^ . (65) 



In the limit of zero temperature, f'(x) — > —5{x) re- 
stricts the k summation to points on the Fermi surface. 
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This allow us to pull one factor of (dE^/de^F = m/m* 
outside the sum. The remaining terms can be computed 
using integration by parts, noting that d 2 ek/dk 2 — 1/m 
and jj J^kn f( E ks _ M) = n s- The resulting expression 
for the conductivity is 



e 
m 



E 



r — iuj{m* /m) s ' 



and when T <^ uj, 



n T (m*/m) i + n^m* /m) ] 



(m*/m)|(m*/m)j 



(66) 



(67) 



The difference 2.M = n-f — nj, between spin up and spin 
down quasiparticlc occupation is obtained by summing 
Eqs. and J2)|). Up to terms of order (B/B ) 3 , 



1 



W 2 



An 2 



f 



B 2 



B 
W 



B 



(68) 



In other words, nj = ^(n+B / Bo) and rif = ^(ti—B/Bq), 
where n = n-\ + = 2 — x is the total Luttinger volume. 
Putting these expressions and Eq. (|62|) into Eq. H67|) gives 



1 1 m ( m* 
— Im- = — 

u> a e n \ m 



1 - (B/B ) 2 



,l + (3/2-2/n)(S/B ) 



(69) 



At large frequencies, the conductivity is related to the 
plasma frequency by 



1 



1 



-— Im — 

lu a 



-in 
7T2"' 



(70) 



Thus, the expression most directly relevant to optical 
measurements is 



4:TTL0 a 



m 
rn 



1 - 



B 



2-i 



(71) 



The coefficient -(5/2 - 2/n) = -[3/2 + x/(2 - x)] is 
negative over the full range of band fillings. 



VI. CONCLUSIONS 

The hybridization mean field theory is appropriate for 
intermediate values of the exchange coupling: to fend off 
competing magnetic states, J must be larger than the en- 
ergy scale associated with spontaneous ordering of the /- 
electron moments; to make contact with realistic metals, 
J should also be less than the conduction bandwidth. In 
this regime, the Kondo and hybridization energies obey 
A K /W ~ V 2 /W 2 ~ a, where a = e ~ 4W ^ 3J is a small 
parameter. The range of valid band fillings, correspond- 
ing to < x < 1 — a, lies between the Kondo- insulator 
and exhaustion limits. 

The hybridization picture captures the essential fea- 
tures of a heavy metal. It allows us to understand how 
quantities such as the effective mass and magnetic sus- 
ceptibility are renormalized by a common enhancement 
factor 1/a. It also explains the existence of the upper 
quasiparticle band and clarifies the role of the hybridiza- 
tion gap. Most important for our purposes, the mean 
field theory provides a simple framework within which to 
compute transport and electrodynamic properties. 

In this paper, we have extended the usual heavy 
fermion picture to the case of an applied Zeeman field. 
The spin splitting of the Fermi surface forces us to dis- 
tinguish between the spin up and spin down quasipar- 
ticles. We have shown that their effective masses differ 
to linear order in B/Bq. Further contributions appear 
at second order because of a reduction in the Kondo en- 
ergy. Accounting for these two effects, we have arrived 
at a specific prediction for the behaviour of the optical 
conductivity. 

Our treatment makes the assumption (jJ,f s < ^2) that 
the applied field is suitably small with respect to the 
Kondo energy: \B\ < 2Akie(1 — x). From an experimen- 
tal point of view, this may still be a large field. The 
Kondo energies of typical heavy fermion materials are 
in the range of 40-100 K. Hence, the predictions in this 
paper may be valid for fields as large as 20-30 T. 

The author would like to thank Dimitri Basov and 
Sasha Dordevic for many helpful discussions. 
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